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ABSTRACT 

We present wide-field niultiband photometry of globular cluster (GC) systems in NGC 1399, NGC 
1404, and NGC 1387 located at the central region of the Fornax galaxy cluster. Observation was 
carried out through J7, B, V , and / bands, which marks one of the widest and deepest f/-band studies 
on cxtragalactic GC systems. The present [/-band photometry enables us to significantly reduce the 
contamination by a factor of two for faint sources (Vq '~23.5). The main results based on some 2000 
GC candidates around NGC 1399, NGC 1404, and NGC 1387 are as follows: (1) the GC system in each 
galaxy exhibits bimodal color distributions in all colors examined, but the shape of color histograms 
varies systematically depending on colors; (2) NGC 1399 shows that the mean colors of both blue 
and red GCs become bluer with increasing galactocentric radius; (3) NGC 1399 shows overabundance 
of GCs in the directions of NGC 1404 and NGC 1387, indicating their ongoing interactions; and (4) 
NGC 1399 also exhibits a ^ 0.'5 offset between the center of the inner GC distribution and the galaxy's 
optical center, suggesting that NGC 1399 is not yet dynamically relaxed and may be undergoing merger 
events. 

Subject headings: galaxies: individual (NGC 1387, NGC 1399, NGC 1404) — galaxies: star clusters: 
general 



1. INTRODUCTION 

The Fornax galaxy cluster is the second nearest clus- 
ter of galaxies after Virgo, with a distanc e modulus of 
(m - M) = 31.51 (jBlakeslee et al.l l2009( ) correspond- 
ing to a distance of 20.0 Mpc. It consists of two ma- 
jor components of the main Fornax cluster centered on 
NGC 1399 and a subcluster dominated by NGC 1316 
(jPrinkwater et an[20bli ). There have been many photo- 
metric and spectroscopic studies of globular cluster (GC) 
systems in the central region of the Fornax cluster, and 
a fair number of them concentrated on the GC system 
associate d with NGC 1 399. 

Bridge&ieFin (|1991[) and IWagner et all (|199H ) pre- 
sented the first CCD photometry for NGC 1399 GCs. 
Both studies presented an unusually high GC specific 
frequency of Sn > 15 (i.e., nu mber of GCs normalized 
to th e galaxy luminosity; see [Harris &: van den Berghl 
19811) Although s ubsequent studies l|Ostrov et al.lll998l: 
Dirsch et al.l 120 03^ derived a lower value {Sn ^ 5), it 



is still highe r than those of normal elliptical galaxies. 
IForbes et afl (jl997f ) suggested that NGC 1404 GCs may 
have been tidally stripped by NGC 1399, resulting in 
a high Sn val ue of NGC 1399 a. n d a l ower Sn value 
of NGC 1404. iKissler-Patig et al.l ^99^ also proposed 
that the over abundance of NGC 1399 GCs can be ex- 
plained by tidal stripping of GCs from neighboring galax- 
ies and by the accretion of G Cs in the gravitationa l 
potential of the Fornax cluster. iGrillmair et al.l (jl999l ) 
supported this idea by showing that GCs in NGC 1399 
and NGC 1404 arc statistically indistinguishable in terms 
of their color distributions and luminosity functions in 



their B, I photometry. iBekki et al.l (j2003D demonstrated 
the tidal stripping and accretion of GCs between NGC 
1399 and NGC 1404 t hrough their numerical simulations. 
iBassino et al.l (j2006al ) presented the asymmetries in the 
azimuthal distribution of NGC 1399 GCs and suggested 
that they may be explained by the tidal stripping of GCs 

from N GC 1387. 

Since lOstrov et all (jl993() first suggested a possible bi- 
modal C—Tl color distribution in the GC system of NGC 
1399, many studies have confirmed the bimodal or multi- 
modal GC color dist ribution in this galaxy through var- 
ious passbands (e.g.. lOstrov et al]ll998l : IGrillmair et al.l 
I1999D . Several studies have in vestigated the prope rties 
of blue and r e d GC s separately. IForbes et al.l (I1998D and 
lOstrov et al.l (|199^ showed that red GCs appea r more 
centrally concentrated than blue GCs. IForte et al.l (j2005[ ) 
suggested that the distribution of blue GCs i s simi lar to 
that inferred for dark matter. iDirsch et al.l (j2003[ ) pre- 
sented that the surface density profile of blue GCs is 
shallower than that of red GCs within 8' but they are 
not distinguishable beyond the radius. They also showed 
that the specific frequency of blue GCs is a fa ctor of three 
larger than that of the red GCs within 7'. IForte et al.l 
(|2005| ) suggested that the specific freque ncies for red and 
blue G Cs are 3.3 and 14.3, respectively. IGrillmair et al.l 
(jl999t) showed that the peak of luminosity function for 
blue G Cs is 0.36 mag bright er in B than that for red 
GCs. iBlakeslee et al.l (j2012D found that the optical- 
infrared color distribution of NGC 1399 GCs is unimodal, 
whereas the optica l colors of the GCs exhibit bimodality. 
IForte et al.l (j2007l ) investigated the connection between 
GCs and the stellar population of NGC 1399 using their 
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TABLE 1 
Host Galaxy Properties 



Galaxy 


R.A. (J2000) 
(h:m:s) 


Decl. {J2000) 
(°:':") 


(mag) 


Type 


NGC 1387 


03:36:56.84 


-35:30:23.85 


12.3 


SBO 


NGC 1389 


03:37:11.67 


-35:44:39.73 


12.8 


SBO 


NGG 1399 


03:38:29.14 


-35:27:02.30 


10.6 


EO 


NGC 1404 


03:38:52.08 


-35:35:37.67 


10.9 


E2 



C — T\ photometry. They suggested that the metal-poor 
subpopulation is homogeneous and exhibits an extended 
spatial distribution contributing ~ 20% of the total stel- 
lar mass, while the metal-rich subpopulation is heteroge- 
neous and dominates the inner region of the galaxy. 

Recent studies based on the Advanced Ca mera for Sur- 
veys (ACS) Fornax Cluster Survey (FCS: [Jordan et all 
|2007[ ) have looked into various properties of GC sys tems 
in the Fornax galaxy cluster. iMieske et al.l ()2010[ ) de- 
tected color-magnitude relation for GCs (i.e., blue tilt) 
and found that the slope is shallower than that found 
for V irgo Cluster Survey GC samples. Master s et alJ 
(poTot ) measured the half-light radii of GCs in ACS FCS 
galaxies, and sugg ested the mean half-light radius as a 
distance indicator. IVillegas et al.l ()2010D found that the 
dispersion of the GC luminosity functions vari es system- 
aticall y with the brightness of parent galaxies. iLiu et aH 
(|2011h investigated the radial color gradients of GC sys- 
tems in early-type galaxies in the Fornax galaxy cluster, 
and found that both red and blue GC subpopulations 
exhibit significant color gradients. They also found that 
the slope of the GC color gradients depends on mass of 
host galaxies. 

There have been a few studies for the GC systems of 
NGC 1404 and NGC 1387. For NGC 1404. lRichtler et alJ 
(|1992f ) presented V , R photometry of the GC system, 
and derived a low specific frequency of Sn = 2.3. They 
suggested that NGC 1404 may have its origin outside the 
Fornax cluster based on the low specific frequency similar 
to that of isolated gal axies together with its high radial 
velocity. iForbes et al.l (fTgOSi) presented broad B ~ I dis- 
tributions for GCs in the inner and outer regions of NGC 
1404 and tentatively suggested a presence of two GC sub- 
populations. The bimodal dist r ibutio n in B — I was later 
confir med by iGrillm air et al.' (1999^ and 'Larse n et alJ 
(|2001h . For NGC 1387, Ba ssino et al . (2006 b) found that 
the C — Tl distribution of the GCs presents a distinct 
separation between red and blue GCs, and the radial dis- 
tribution of red GCs is more centrally concentrated than 
blue GCs. 

Several studies have been done on kinematics and 
dynam ics for GCs around NGC 1399. IGrillmair et al.l 
(|1994f l suggested that, based on the low-dispersion spec- 
tra of 47 GCs, the large velocity dispersion of the GCs 
indicates that they are associated with the gravitational 
potential of t he Fornax cluster rather t han with that of 
NGC 1399. iKissler-Patig et al.l (IT999I ) supported this 
suggestion by finding a strong radial increase of the 
velocity dispersion for 74 GCs. They found no com- 
pelling evidence for rotation and could not detect differ- 
ences betweeii the k inematics of the blue and red GCs. 
iRichtler et al.l ()2004[ ) obtained medium-resolution spec- 
tra of 468 GCs, and they did not find any signature of 



TABLE 2 
Observation Log 



Date (UT) 


Filter 


Exposure Time 


Number of 


Seeing 






(s) 


Frames 


(arcsec) 


2006 Nov 24 




1800 


6 


0.81-1.08 


2006 Nov 24 


B 


900 


3 


0.91-1.20 


2006 Nov 24 


V 


240 


3 


0.94-1.10 


2006 Nov 24 


I 


300 


6 


0.63-0.80 


2006 Nov 28 


[/a 


1800 


9 


1.14-1.43 


2006 Nov 28 


B 


900 


3 


1.15-1.35 


2006 Nov 28 


V 


360 


3 


0.86-1.07 



^ We used the SDSS u filter and later ealibrated the photometry 
to the Johnson U system. 

rotation. They suggested that the large velocity disper- 
sion of the GCs co uld arise from th e cont amination by 
foreground objects. iSchuberth et al.l ()2010[ ) investigated 
kinematical and dynamical properties of the GC system 
of NGC 1399 using the spectra of some 700 GCs. They 
found that blue and red populations are kinematically 
distinct, and the properties of red population resemble 
those of field stellar populations. 

Several spectroscopic studies investigated chem- 
ical abundance a nd ag es of NGC 1399 GCs. 
iKissler-Patig et ahl ()1998l) obtained spectra of 18 
GCs in NGC 1399 and found that their abundances 
span the range observed in the Milky Way and M31. 
They pointed out that V — I correlates well with 
metallicity but the slope of the relation becomes flatter 
toward redder colors if the values of Milky Way GCs are 
considered together. iForbes et al.l (j2001f ) investigated 
the ages of the GCs and found that the majority of their 
sample (10 GCs) have ages of around 11 Gyr. 

In this paper we provide a new [/, B, V, and / pho- 
tometric catalog of GC candidates in the central region 
of the Fornax galaxy cluster. This is one of the widest 
and deepest J7-band studies on extragalactic GC sys- 
tems. Section 2 describes the observations and data 
reduction. Section 3 gives the GC catalog for NGC 
1399, NGC 1404, and NGC 1387, discussing the selection 
method and the foreground and background contamina- 
tion level. Section 4 presents the properties of GC sys- 
tems around the galaxies, including spatial distributions, 
color-magnitude diagrams (CMDs), color distributions, 
and the radial variations of color bimodality. Section 5 
summarizes our results. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observations 

The observations were carried out on the nights of 2006 
November 24 and 28 (UT) with the Mosaic II CCD im- 
ager mounted on the prime focus of the 4 m Blanco 
telescope at Cerro Tololo Inter-American Observatory 
(CTIO). The observing conditions were clear and photo- 
metric with an average seeing of around 1". The Mosaic 
II consists of eight 2048x4096 pixel CCDs with a pixel 
scale of 0'.'27 providing a total field of view of 36'x36'. 

Our target field was chosen such that the brightest 
elliptical galaxy in the Fornax cluster, NGC 1399, and 
its three neighboring galaxies, NGC 1404, NGC 1387, 
and NGC 1389, were all imaged in a single telescope 
pointing. Figure [T] shows the median combined image of 
our target field as imaged in the i?-filter, and Table[T]lists 
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Chip 4 



Chip 3 



NGC^l-387 



NGC 1404 



Chip 6 



Chip 2 



NGC 1.389 



Chip 5- 



Chip 1- 



Fig. 1. — Median combined S-band image of our target field with tlie configuration of Mosaic II chips overlaid. The gaps between chips 
were filled in by creating mosaics using a three-point dither pattern. The image covers an area of 36' X 36', containing four Fornax cluster 
galaxies as labeled and their GC systems. North is up, and east is to the left. 



the basic properties of the four members of the Fornax 
cluster. We observed the target field with u, B, V, and 
/ filters. We note that the standard Johnson U filter 
was not available at the time of our observations, so we 
used the Sloan Digital Sky Survey (SDSS) u filter as an 
alternative and later calibrated the photometry to the 
Johnson U system. The resulting effect of this calibration 
procedure is discussed below. The total exposure times 
in u, B, V, and / were 27000, 5400, 1800, and 1800 s, 
respectively. We split our observations in each band into 
two sets of three dithered exposures (five sets for u) to 
fill in the gaps between the individual CCD chips in the 
mosaic, and to minimize the impact on our photometry 
from CCD blemishes (e.g., hot pixels and bad columns). 



A summary of observations is provided in Table [2] 

2.2. Pre-processing and Photometry 

All im ages were pr e-processed using the MSCRED 
package (iValdesI [19981) in IRAF0 We used the CCD- 
PROC task to correct for cross-talk between the CCD 
chips, apply and trim the overscan, correct for the 
bias level, and apply the flat-field correction. We used 
MSCPIXAREA to correct for the variations in pixel 
scales across the frame produced by geometric distortion 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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o 2006 Nov. 
X 2007 Mar. 




-2-10 1 2 3 

U-B 

Fig. 2. — Difference between the standard and instrumental U magnitudes of standard stars after photometric zero point and atmospheric 
extinction corrections, as a function of U — B. In the calibration of U band, all the standard stars observed in subsequent observing 
runs during three years were used to derive a reliable solution. Different symbols are used to distinguish observations from different 
seasons. Two dashed lines represent error-weighted least-squares fits for blue and red standard stars, respectively. Red solid line repre- 
sents the adopted best-fit function, which connects the two functions using internal division points of them between two vertical orange lines. 



in the Mosaic II imager. Each mosaic image was then 
split into eight individual images using the MSCSPLIT 
task. 

To produce a clean high signal-to-noise image in each 
filter, we combined the images usin g the M0NTA GE2 
routine in the ALLFRAME packag e (IStetsonlll994[) . We 
then applied a ring median filter (|Seckerill995T having 
an inner radius of 2.5 times the mean full width at 
half-maximum (FWHM) and a width of 5 pixel on the 
combined images and subtracted off these images from 
the original images to remove the rapidly varying back- 
ground light of the galaxies while preserving the point- 
sources. Objects were detected from the background- 
subtracted images using the DAOPHOT II/ALLSTAR 
routines through the following steps. First, we derived 
a list of objects and created new images with these ob- 
jects subtracted. From these images, we detected fainter 
objects that were missed in the previous step and cre- 
ated another set of object-subtracted images. This pro- 
cedure was iterated five times, and all detected objects 
were merged into a master catalog. The master catalog 
was then input to the ALLFRAME program (jStetsoni 
|1994[ ) along with all of the images and their point- 
spread function (PSF) models for final PSF photome- 
try. Aperture correc tions were co mputed through the 
DAOGROW routine (|Stetsonlll990h using bright and iso- 



lated point-sources, and we applied them to our photom- 
etry above. Finally, an error- weighted mean instrumen- 
tal magnitude for each object was obt ained using t he 
DAOMATCH/DAOM ASTER routines (|Stetsonlll993l ). 

The photometric calibration was achieved using stan- 
dard stars in the fields of [Landolt ( 1992) and Stetso^ 
([20M l. During the course of our observing runs, we im- 
aged eight standard fields placing the stars of interest 
on chip 6. Each image contains 15-25 standard stars, 
which covers an appropriately large range in color. We 
also observed the standard fields in many different air- 
masses. In the calibration of U band, however, we used 
all the standard stars that are observed in subsequent 
observing runs during 2006, 2007, and 2008, because the 
photometric uncertainty in U band is quite large and the 
number of standard stars is too small to derive reliable 
solutions. To calibrate the instrumental magnitudes to 
the Johnson-Cousins standard magnitudes, we used the 
transformation equation of the form 

ACtd = Wi„rt +cl -f c2 X airmass + c3 x color + c4 x color^, 

where M^^^ is the calibrated magnitude, rriinst is the in- 
strumental magnitude that is normalized to 1 s, cl is the 
photometric zero point offset, c2 is the extinction coef- 
ficient for airmass corrections, and both c3 and c4 are 
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TABLE 3 

Photometric Calibration Coefficients 



Date (UT) Filter cl c2 c3 c4 



2006 Nov 24 C/biuc=' -2.1194 ± 0.0127 -0.4495 ± 0.0072 0.1458 ± 0.0038 NaN 

C/,od=' -2.0699 ± 0.0380 -0.4442 ± 0.0144 -0.0819 ± 0.0581 0.1390 ± 0.0259 

B 0.3711 ± 0.0045 -0.2236 ± 0.0037 0.0803 ± 0.0015 NaN 

V 0.6402 ± 0.0035 -0.1386 ± 0.0028 -0.0297 ± 0.0011 NaN 

/ 0.0578 ± 0.0039 -0.0576 ± 0.0030 -0.0099 ± 0.0011 NaN 

2006 Nov 28 C/biuo'' -2.1641 ± 0.0128 -0.4495 ± 0.0072 0.1458 ± 0.0038 NaN 

(7,<3d^ -2.1196 ± 0.0383 -0.4442 ± 0.0144 -0.0819 ± 0.0581 0.1390 ± 0.0259 

B 0.3243 ± 0.0048 -0.2236 ± 0.0037 0.0803 ± 0.0015 NaN 

y 0.6061 ± 0.0037 -0.1386 ± 0.0028 -0.0297 ± 0.0011 NaN 

/ 0.0350 ± 0.0041 -0.0576 ± 0.0030 -0.0099 ± 0.0011 NaN 



For (7 band, two functions are used to find the best- fit (see the text). 



TABLE 4 

Chip-to-chip Correction Coefficients 



Chip 


Filter 


cl 


c2 


Chip 1 


U 


0.0347 ± 0.0020 


-0.0334 ± 0.0055 




B 


0.0236 ± 0.0007 


-0.0147 ± 0.0022 




V 


-0.0165 ± 0.0004 


0.0029 ± 0.0014 




I 


0.0037 ± 0.0005 


0.0010 ± 0.0017 


Chip 2 


u 


0.0323 ± 0.0019 


-0.0609 ± 0.0055 




B 


0.0162 ± 0.0007 


-0.0094 ± 0.0024 




V 


0.0034 ± 0.0005 


0.0046 ± 0.0015 




I 


0.0023 ± 0.0005 


-0.0037 ± 0.0017 


Chip 3 


u 


0.0087 ± 0.0019 


-0.0675 ± 0.0057 




B 


0.0197 ± 0.0007 


-0.0125 ± 0.0021 




V 


0.0020 ± 0.0004 


0.0038 ± 0.0014 




I 


0.0085 ± 0.0005 


0.0013 ± 0.0016 


Chip 4 


u 


0.0208 ± 0.0020 


-0.0156 ± 0.0062 




B 


0.0195 ± 0.0007 


-0.0093 ± 0.0021 




V 


-0.0060 ± 0.0005 


0.0016 ± 0.0015 




I 


-0.0055 ± 0.0005 


0.0026 ± 0.0016 


Chip 5 


u 


0.0409 ± 0.0020 


0.0013 ± 0.0057 




B 


-0.0024 ± 0.0007 


-0.0018 ± 0.0022 




V 


-0.0149 ± 0.0005 


0.0022 ± 0.0016 




I 


0.0019 ± 0.0007 


-0.0006 ± 0.0014 


Chip 7 


u 


-0.0183 ± 0.0020 


-0.0016 ± 0.0060 




B 


-0.0052 ± 0.0007 


0.0047 ± 0.0021 




V 


-0.0023 ± 0.0004 


-0.0002 ± 0.0014 




I 


0.0031 ± 0.0005 


0.0033 ± 0.0017 


Chip 8 


u 


0.0136 ± 0.0020 


-0.0069 ± 0.0060 




B 


0.0083 ± 0.0007 


-0.0037 ± 0.0021 




V 


-0.0046 ± 0.0005 


0.0005 ± 0.0015 




I 


0.0122 ± 0.0005 


0.0035 ± 0.0017 



the color coefficients. The colors used in the calibration 
axeU - B, B - V, B - V, &nAV - I for U, B, V, and 
/, respectively. The photometric coefficients are deter- 
mined by an iterative process of error-weighted fitting of 
the transformation equation to the standard magnitudes. 
The quadratic color term is included only in the equation 
of U band, because a significant deviation from a linear 
fit occurs in the transformation of U band. Figure [2] 
shows the U magnitude difference between the standard 
and instrumental magnitudes corrected for both of the 
photometric zero point and the atmospheric extinction, 
as a function oi U — B. The dashed lines represent two 
best-fit functions for blue and red standard stars, respec- 



tively. However, if we divide the observed data into two 
groups by U — B and use the two functions in calibration 
separately, the discontinuity between two functions may 
lead to an artificial clump in color distributions of the 
data. To avoid this problem, we connected the two func- 
tions (red lines) at the internal division points of them 
in the range 0.29 < U — B < 1.35 (two vertical lines). 
Finally, the instrumental magnitudes were converted to 
the standard Johnson-Cousins UBVI system using the 
set of photometric coefficients listed in Table |3] 

Calibrating the SDSS u photometry to the Johnson U 
system may cause systematic shifts in the colors based 
on this band since the system throughputs arc different 
between the two filters. To estimate these shifts, we sim- 
ulated our observations using the system throughputs of 
the two filters, synthetic sp ectra of GCs from the YEPS 
model (jChung et al.l [20121 ). and stars from the Kurucz 
ATLAS9 model. We calculated the magnitude through 
each filter for each model GC, and converted it to the 
standardized Johnson U magnitude using the synthetic 
spectrum of a star having the same color of the GC. We 
then calculated the differences between these standard- 
ized magnitudes from the respective filters. The differ- 
ences range from 0.002 mag to 0.022 mag, depending on 
the colors of the model GCs. We conclude that these are 
small enough to have negligible impact on our photomet- 
ric analyses throughout the paper. 

As the individual CCD chips of the Mosaic may have 
different color terms, we attempted to correct for this 
"chip-to-chip variation" as follows. On the night of 2008 
May 5, we took multiple uBVI observations of the Lan- 
doh standard fields SA98 and SA107 with the Mosaic 
II, such that the same set of standard stars was posi- 
tioned on each chip. Using these data, we derived the 
color terms that correct for the chip-to-chip variation. 
Although these color terms were very small, we applied 
them to our data in hopes of making the multiple mea- 
surements of the same source more consistent. The equa- 
tion used in correction is of the form 

Af,„, = m + cl + c2 X color, (2) 

where M^„^ is the magnitude converted to that on the 
reference chip 6, m is the magnitude measured on each 
chip, cl is the zero point offset, and c2 is the color co- 
efficient. The colors used in the calibration arc U — B, 
B~V, B~V, and y — / for U, B, V, and /, respectively. 
Table [4] lists the zero point offsets and color coefficients 
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F ig. 3. — Completenes s functions obtained from artificial star tests. Detection efficiencies (filled circle) were fitted to the analytic function 
by IFlemine et aU I I1995I) . which arc denoted by solid lines. The purple, blue, green, and red lines represent the U, B, V, and / bands, 
respectively. The functions are for the overall image, except the galactic central regions and the gaps between mosaic chips where the 
detection efficiency rapidly drops. Limiting magnitudes where the completeness drops to 50% are indicated in the upper-right corner. 



for each chip. 

We corrected for the foreground Gala ctic extinction 
using the reddening maps provided by iSchlegel et al.l 
(|1998f ) . To account for the atmospheric dispersion correc- 
tor (ADC) of the 4 m Blanco telescope that significantly 
modifies the system throughputs, we computed the effec- 
tive wavelengths of our four bands from the total system 
throughputs including the effect of the ADC, and we de- 
rived A\/E {B — V) values based on the extinction law 
provided bv lCardeUi et all (fT989t) . The mean extinction 
values in J7, B, V, and I are 0.06, 0.05, 0.04, and 0.02 
mag, respectively. 

Our images were calibrat ed for astrometry using the 
USNO-B 1.0 catalog stars (IMonet et al.ll2003D and the 
MSCTPEAK task. The typical rms of our astromctric 
solutions is at the level of 0'.'2. The R.A./decl. coordi- 
nates of our subsequent source catalogs use the positions 
derived from these images. 

2.3. Completeness 

We performed artificial star tests to estimate the com- 
pleteness of our photometry. We used the ADDSTAR 
routine of DAOPHOT II to generate 4000 artificial stars 
in each chip with magnitudes and color ranges consis- 
tent with those of the point-sources detected and mea- 
sured above. The luminosity function of our artificial 
stars was chosen to be flat to preserve statistical signifi- 
cance in all magnitude ranges. Artificial stars were added 



to the PSF-subtractcd images using the empirical PSFs 
derived during our photometry. We used the same pho- 
tometry procedure as above to detect and measure arti- 
ficial stars in these images. This process was repeated 
seven times with the positions of artificial stars altering 
randomly each time, resulting in a total of 28,000 arti- 
ficial stars per chip. In each 0.3 mag bin, we calculated 
the detection efficiency, i.e., the ratio of the number of 
recovered artificial stars to the number of added ones. 
Central regions of galaxies and the gaps between CCD 
chips were masked for the calculation. To obtain a limit- 
ing magnitude at which the completeness drops to 50 %, 
the det ection efficienci e s wer e fitted to the analytic func- 
tion bv [Fleming et al.l (II 9951 ). 

/M-^(i- (3) 

y 1 + a (to — 7Ti,imj^ 

where m,i„ is the limiting magnitude, and a is the pa- 
rameter that controls the steepness of the completeness 
function near the limiting magnitude. Figure [3] shows the 
measured completeness (dots) and the fitted functions 
(solid lines) in U, B, V, and / bands depicted by purple, 
blue, green, and red colors, respectively. The limiting 
magnitudes of the four bands are 24.4, 24.4, 23.7, and 
22.8 mag in U, B, V, and / bands, respectively. We note 
that the limiting magnitudes slightly vary in the < 0.2 
mag level among individual CCDs due to variations in 
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Fig. 4. — y-band completeness functions at various radial bins from the center of galaxies — NGC 1399 (top), NGC 1404 (middle), and 
NGC 1387 (bottom). The black thick line represents the overall completeness function. The completeness is very low near the center of 
the galaxies, and increases as one moves out from the center, and then flattens out beyond 1.'3 for NGC 1399, O'.l for NGC 1404, and 0'.5 
for NGG 1387. 




Fig. 5. — Photometric error and image parameters (x and Sharpness) in the U band as a function of magnitude. Sohd hnes show the 
point-source selection criteria. Black dots represent the selected sources with good photometry, while gray dots represent the sources that 
do not satisfy any of the selection criteria. 
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detector characteristics. 

We also investigated the radial variations of complete- 
ness around NGC 1399, NGC 1404, and NGC 1387. In 
the central region of a galaxy, the detection effectiveness 
at a certain magnitude decreases as the surface bright- 
ness of the galaxy increases, and it finally becomes zero 
at a saturated region. To quantify this effect in our obser- 
vations, we performed separate artificial star tests near 
the galaxy centers. We added 8000 artificial stars on the 
individual CCD chips that contain the galaxies (NGC 
1399, NGC 1404, and NGC 1387), this time with a cen- 
trally concentrated radial distribution centered on each 
galaxy. This experiment was repeated five times with the 
positions of artificial stars altering each time. The com- 
pleteness was then calculated as a function of magnitude 
in concentric circles with radial intervals of 0'.2. Figure |4] 
shows the V-band completeness functions at various ra- 
dial bins from the center of each galaxy. The thick black 
line represents the completeness function for the overall 
region, as derived above. As expected from the strong 
background lights near the centers of our target galaxies, 
we find that the completeness is very low near the cen- 
ter. The completeness increases as one moves out from 
the center, and then flattens out beyond 1^3 for NGC 
1399, 0.'7 for NGC 1404, and 0'.5 for NGC 1387. 

3. CONSTRUCTION OF GLOBULAR CLUSTER CATALOGS 
3.1. Selection of Globular Cluster Candidates 

At the distance of the Fornax cluster 20.0 Mpc; 
iBlakeslee et "all [20091) . 1 pc is equal to O'.'Ol, so a typi- 
cal GC appears as a point-source in our Mosaic images. 
We therefore can not select GCs based on their apparent 
sizes but must rely on photometric colors. Details on the 
selection process of GC candidates are given below. 

We start by using the ALLFRAME parameters CHI 
and SHARP to select point-sources from our photometric 
catalog. Figure [5] shows the SHARP, CHI, and the ALL- 
FRAME photometric error as a function of Uq magni- 
tude with solid lines showing the selection criteria. While 
the figure only shows the [/-band selection, we applied 
the same selection criteria to all bands simultaneously, 
and the resulting selection includes 5500 point-source ob- 
jects. We then employ color cuts using our multi band 
photometry as demonstrated in Figure |6] to select GC 
candidates. The blue selection boxes shown in the three 
color-color diagrams were determined on the basis of the 
findings below. First, we matched our point-source ob- 
j ects w ith the spectroscopic catalog of iSchuberth et ahl 
((2OT0h and found that the NGC 1399 GCs, as identified 
by their measured line-of-sight velocities, form a tight se- 
quence (red dots in Figurc|6l) in the color-color diagrams, 
whereas foreground stars are located along an extended 
sequence (cyan dots in Figure |6]). Second, we downloaded 
Hubble Space Telescope {H ST) /ACS images of NGC 
1399 from the archive and identified genuine GCs based 
on their sizes . The sizes were measured using ISHAPE 
(jLarse h 1999") and fo ll owing the same procedure out- 
lined in Strader et al.l ()2006[ ). We then cross- identified 
the GCs with our point-source objects, and found that 
most of them lie inside the tight sequence (orange dots 
in Figure in])- To select GC candidates, we further ap- 
plied a magnitude cut of Vb < 20.6 (corresponding to 
My < —11) to the color-color-selected sample. This 




0.0 0.5 LO 1.5 2.0 

(V-I)o 




0.0 0.5 1,0 1.5 3.0 

(V-I)o 




1.5 



1 2 3 

(B-I)o 

Fig. 6.— Color selection of GC candidates in NGC 1399. The 
selection criteria demarcated by blue lines were determined based 
on confirmed GCs wi th radial velocity measurements (red dots; 
ISchuberth et all 120101 ) and/or with size measurements (orange 
dots; derived from HST ACS data). The cyan asterisks repre- 
sent foregro und stars confirmed with radial velocity measurements 
ISch uberth e t al.ll201Cll l. GC candidates are clearly separated from 
background sources in the {U — B)q color. 
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magnitude limit was determined base d on the bright- 
est G C in the spectroscopic catalog of iSchuberth et ahl 
()2010[) . and we note that it is consistent with the hmit 
used bv lMieske et al.l ([2006 ) to separate GCs from their 
ultra compact dwarf galaxies sample. The color-color 
and magnitude cuts result in a final sample of 2037 GC 
candidates. 

Figure |6] demonstrates the advantage of using U band 
for separating star-forming background galaxies from the 
extragalactic GCs. The background galaxies have {B — 
y)o colors similar to the bluer GCs, but due to their star- 
forming nature, galaxies exhibit significantly bluer colors 
in {U — V)q and {U — B)q than the GCs. This allows us to 
effectively reduce contaminati on by background ga laxies. 
The same effect was found bv lDirsch etall (I2003D using 
the C — n Washington color index. Nonetheless, our 
selection box still contains several foreground stars, and 
it is important to quantify the contamination level of our 
color-color-selected objects. We provide the estimation 
of the contamination level below. 

3.2. Estimation of Contamination 

To estimate the number of contaminating sources (fore- 
ground stars and background galaxies) in our GC can- 
didate list, it would be best to use a control field ob- 
served near the science field with the same observational 
conditions. However, this approach is impractical for us 
because the J7-band observations require very long expo- 
sure times. As an alternative, we used two independent 
methods to estimate the contamination level of our GC 
candidate list. 

On average, chip 1 is located farthest away from the 
target galaxies NGC 1399, NGC 1387, and NGC 1404, 
so it has the least number of GCs expected among the 
eight CCD chips In Figure [3 we compare the {U — 
B)o versus {B — /)o color-color diagrams of the 5500 
point-source objects for each chip. In chip 1, there are 
41 objects that are selected as GC candidates using the 
color-color selection box as in Figure [51 If we assume 
that all 41 objects are contaminants, this implies that 
our catalog of GC candidates has a contamination level of 
16 %. We note that 16 % is in fact an upper limit because 
chip 1, which is at an average distance of 24' 140 kpc) 
from the center of NGC 1399, is likely to contain several 
GCs given that the NGC 1399 GCs are known to extend 
out to a projected d istance of at least 45' 260 kpc; 
iBassino et alll2006aD . 

We further estimated the number of foreground stars 
using the B esancon model of the Milky Way population 
([Robin et al. 2003). We generated a synthetic catalog 
of stars in the direction of the sky consistent with our 
target field and then selected stars using the same color- 
color cut as used in Figure [S] The estimated number of 
foreground stars in our GC candidate list turned out to 
be 180 and the corresponding contamination level is 9 %. 
As expected, this is lower than the upper limit set by 
using chip 1 as the control field above. Based on the two 
estimates, the contamination level of our GC candidate 
list is likely to be in the 9 % -16 % range. 

Although the contamination level itself may not be 
that high, it is important to test whether the contaminat- 

^ NGC 1389 is known to have very few GCs (see Section 4.1), so 
we only considered the other three Fornax galaxies here. 



ing sources will affect our GC color distribution analyses. 
Figure [H shows the comparison of {U — B)o color dis- 
tributions for all point-sources (shaded histogram), our 
GC candidate list (blue), and synthetic stars from the 
Besangon model (red) between chips 7 and 1. Based on 
the similarities found in the red and blue color distribu- 
tions in the upper panel of Figure [H chip 1 seems to be a 
good control field. From the same figure, we find that the 
foreground MW stars are expected to have a. {U — B)q 
color distributions that is peaked toward the bluer side 
similar to the {U — B)q color distribution of GC can- 
didates. However, as demonstrated in the lower panel 
of Figure [51 the number of these foreground stars, when 
compared to that of the GC candidates, is too small to 
have any impact on our color distribution analyses pre- 
sented in Sections 4.3 and 4.4. 

3.3. Photometric Catalog 

Table [5l presents the photometric catalog of total 2.037 
GC candidates. The columns give the following informa- 
tion. 

Column 1. Identification number 

Columns 2 and 3. Right ascension and declination 
(J2000) 

Columns 4 and 5. Uq magnitude and its rms uncertainty 
Columns 6 and 7. Bq magnitude and its rms uncertainty 
Columns 8 and 9. Vq magnitude and its rms uncertainty 
Columns 10 and 11. Iq magnitude and its rms uncer- 
tainty. 

We note that all magnitudes in the catalog arc cor- 
rected for Galactic extinction. 

4. RESULTS AND DISCUSSION 
4.1. Spatial Distribution 

Figure [9l maps the spatial distributions of GC can- 
didates (left panel) and foreground/background objects 
(right). The color-filled contours show the distribution of 
surface number density. The black crosses mark the po- 
sitions of the optical centers of galaxies (from top, NGC 
1399, NGC 1387, NGC 1404, and NGC 1389). In the left 
panel, three GC systems are clearly visible around NGC 
1399, NGC 1404, and NGC 1387. The number of GCs 
around NGC 1389 is too small, so we omit the galaxy 
from the following analysis. In the right panel, the distri- 
bution of foreground/background objects appears fairly 
uniform. The slight overdensity in the region between 
NGC 1399 and NGC 1387 is likely due to the presence of 
a distant galaxy cluster, given that many objects in this 
region are bluer than -0.3 in {U — B)q. Another over- 
dense region is right on the center of NGC 1399, which 
indicates that GCs in highly dense regions can be mis- 
classified as foreground/background objects in our GC 
selection procedure, as mentioned in Section 3. Note that 
the residuals of the mosaic CCD gaps are discernible as 
a vertical feature in the middle of the panel. 

It is noteworthy that the GC distribution around NGC 
1399 shows interesting features. On the one hand, there 
is an overabundance of GCs in the outer region of NGC 
1399 toward N GC 1404 and NGC 13 87. This was also 
pointed out bv IBassino et al.l ([2006af ). and it may indi- 
cate interactions of NGC 1399 with NGC 14 ( 34 and NGC 
1387 in the recent past ([Forbes et al.lll997t iBekki et all 
|2003() . On the other hand, a ~ 0^5 displacement is found 
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Fig. 7. — Color-color diagrams of point-sources on each chip. Solid boxes indicate the selection criteria for the GC candidates and the chip 
numbers arc indicated in the lower-left corner of each panel. Chips 7, 6, and 3 contain NGC 1399, NGC 1404, and NGC 1387, respectively. 
Chip 1, which is located farthest away from the galaxies, is used to estimate the upper limit of the contamination level of our data. 



between NGC 1399's optical center and the center of the 
GC distribution in the inner region of the galaxy as de- 
noted by a red contour in the density map. iPaolillo et al.l 
(|2002f ) reported an asymmetric X-ray halo for NGC 1399, 
consisting of three components with different centers; 
a central component, a galactic component centered 1' 
southwest of NGC 1399, and a cluster component cen- 
tered ~5f6 northeast of the galaxy. While the center 
of the GC distribution does not coincide with the X-ray 
centers, the asymmetric distribution of the GCs may sug- 



gest that the NGC 1399 is not yet dynamically relaxed 
and may be undergoing merger events. 

Figurc[TO]presents the surface number densities against 
the galactoccntric radius for the three galaxies. In order 
to examine the radial extent of each GC system, we cal- 
culated the surface GC number density as a function of 
the distance from the galaxy center. We first set up a 
series of radial bins (annulus), each containing approxi- 
mately equal numbers of GCs (- 100 GCs for NGC 1399; 
-40 GCs for NGC 1404; -25 GCs for NGC 1387). The 
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Fig. 8. — U — B color distributions for all point-sources (shaded histogram), GC candidates (blue), and synthetic stars from Besangon 
model (red) in chip 1 and chip 7. The similarities found in the red and blue color distributions in the upper panel indicate that chip 1 is a 
good control field. 
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Fig. 9. — Spatial distribution of GC candidates (left) and foreground/background objects (right) over color-filled density contour maps 
with the color scale from black to red for increasing surface number density. Black crosses represent the positions of the optical centers of 
galaxies (from top, NGC 1399, NGC 1387, NGC 1404, and NGC 1389). 



TABLE 5 

UBVI Photometry of GC Candidates in the Fornax Galaxy Cluster 



GC ID 


R.A. (J2000) 


Decl. (J2000) 


Uo 


U Error 


Bo 


B Error 


Vo 


V Error 


lo 


/ Error 


1 


54.4670781 


-35.8412594 


23.863 


0.057 


23.466 


0.018 


22.611 


0.025 


21.614 


0.018 


2 


54.4048663 


-35.8368808 


23.099 


0.028 


22.941 


0.011 


22.069 


0.015 


21.119 


0.010 


3 


54.1532595 


-35.8313662 


23.048 


0.038 


22.743 


0.012 


21.970 


0.015 


20.920 


0.020 


4 


54.3887923 


-35.8307234 


24.085 


0.070 


23.983 


0.029 


23.384 


0.049 


22.462 


0.030 


5 


54.3871824 


-35.8300091 


23.772 


0.054 


23.865 


0.025 


23.147 


0.038 


22.327 


0.028 


6 


54.3545131 


-35.8284302 


23.147 


0.030 


22.633 


0.008 


21.661 


0.011 


20.490 


0.006 


7 


54.5147090 


-35.8246663 


22.822 


0.023 


22.949 


0.014 


22.247 


0.015 


21.430 


0.013 


8 


54.3057038 


-35.8223111 


24.007 


0.061 


24.027 


0.028 


23.363 


0.045 


22.429 


0.028 


9 


54.2008513 


-35.8137756 


24.064 


0.070 


24.111 


0.030 


23.383 


0.045 


22.486 


0.030 


10 


54.4062928 


-35.8090401 


21.291 


0.007 


21.338 


0.004 


20.674 


0.005 


19.854 


0.003 



Note. The first 10 rows of the table arc reproduced here; the full version of this table containing the 2,037 GC candidates is 
available online at the CDS. 
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Fig. 10. — Number of GC candidates per unit area projected on the sky as a function of radial distance from the center of the galaxy. 
Radial bins are defined to contain approximately equal numbers of GCs (~ 100 GCs for NGC 1399; ~ 40 GCs for NGC 1404; ~ 25 GCs 
for NGC 1387). The numbers are corrected for the spatial completeness, but not for the contamination by foreground/background objects 
(see the text). Dashed line shows de Vaucouleurs' i? fit of each profile. 



number of GCs in each bin was then corrected using the 
corresponding completeness function of the bin. We fi- 
naUy divided the corrected numbers of GCs by the area of 
annuli to get surface densities. The surface density pro- 
file of each galaxy is well described by de Vaucouleurs' 
i?^/* fits (dashed lines). We determine the limiting radii 
of GC systems of NGC 1404 (3:4) and NGC 1387 (3.7) 
as the radial points where the density profiles begin to 
depart from de Vaucouleurs' law. Since the GC systern . 
of NGC 1399 is extended (-45'; iBassino et all l2006aD 
well beyond the field of view (36' x 36') of our Mosaic 
observations, we regard the entire GCs in the images as 
NGC 1399's GC system, except for the regions inside the 
limiting radn of NGC 1404 and NGC 1387. 



Figure [TT] presents the CMDs for all point-sources. 
Magnitudes and colors were corrected for the Galac- 
tic extinction, as mentioned in Section 2.3. The black 
and gray dots represent the GC candidates and fore- 
ground/background objects, respectively. The dashed 
lines indicate the limiting magnitude of 50 % complete- 
ness. The tilt of the limiting magnitude lines results from 
lower detection efficiency of red GCs than blue GCs in 
the U band. The {U — B)o colors corresponding to GCs, 
background galaxies, and foreground stars can be deter- 
mined both by the {U — B)o versus {B — I)o diagram 
(the bottom panel of Figure [6|) and by the Vq-{U — B)o 
CMD (the left panel of Figure [TT]) . GCs are placed at 
—0.2 < {U — B)o < 0.8, while most background galaxies 
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Fig. 11. — Color-magnitude diagrams of point-sources. Black dots represent GC candidates, and gray dots represent fore- 
ground/background objects. Observational uncertainties as a function of Vb are shown by error bars in the left corner of each diagram. 
Dashed lines indicate the limiting magnitude. The tilt of the limiting magnitude results from lower detection efficiency of red GCs than 
blue GCs in the U band. The GC candidates fainter than the limiting magnitude are taken out from the subsequent analysis. 



TABLE 6 

Results from GMM Analysis for Color Distributions of GCs 



Color 






fir 




fr 


ar/o-f, 




p{DD) 


p{kurt) 


NGC 1399 


U 


- B 


0.016 ± 0.006 


0.102 ± 0.003 


0.372 ± 0.021 


0.188 ± 0.009 


0.34 


1.84 


0.010 


0.090 


0.860 


U 


- V 


0.670 ± 0.007 


0.113 ± 0.005 


1.147 ± 0.022 


0.269 ± 0.008 


0.42 


2.38 


0.010 


0.100 


0.140 


u 


- I 


1.557 ± 0.008 


0.134 ± 0.006 


2.188 ± 0.025 


0.357 ± 0.008 


0.45 


2.66 


0.010 


0.080 


0.010 


B 


- V 


0.673 ± 0.004 


0.058 ± 0.002 


0.849 ± 0.007 


0.060 ± 0.003 


0.34 


1.03 


0.010 


0.070 


0.010 


B 


- 1 


1.571 ± 0.005 


0.087 ± 0.003 


1.925 ± 0.014 


0.135 ± 0.007 


0.37 


1.55 


0.010 


0.050 


0.010 


V 


- 1 


0.910 ± 0.005 


0.063 ± 0.003 


1.108 ± 0.013 


0.068 ± 0.007 


0.30 


1.08 


0.010 


0.060 


0.010 



NGC1404 



u - 


B 


-0.004 ± 0.012 


0.079 ± 0.008 


0.490 ± 0.024 


0.120 ± 0.018 


0.39 


1.52 


0.001 


0.017 


0.001 


u - 


V 


0.682 ± 0.015 


0.101 ± 0.011 


1.393 ± 0.026 


0.136 ± 0.017 


0.39 


1.35 


0.001 


0.001 


0.001 


u - 


I 


1.631 ± 0.018 


0.124 ± 0.015 


2.583 ± 0.029 


0.161 ± 0.019 


0.39 


1.30 


0.001 


0.001 


0.001 


B - 


V 


0.686 ± 0.007 


0.052 ± 0.007 


0.901 ± 0.005 


0.029 ± 0.004 


0.38 


0.56 


0.001 


0.014 


0.001 


B - 


I 


1.634 ± 0.011 


0.075 ± 0.011 


2.091 ± 0.009 


0.051 ± 0.006 


0.39 


0.68 


0.001 


0.001 


0.001 


V - 


I 


0.949 ± 0.006 


0.041 ± 0.006 


1.190 ± 0.006 


0.030 ± 0.005 


0.39 


0.73 


0.001 


0.001 


0.001 



NGC1387 



u 


~ B 


-0.019 ± 0.017 


0.084 ± 0.011 


0.446 ± 0.030 


0.128 ± 0.017 


0.49 


1.52 


0.001 


0.041 


0.001 


u 


- V 


0.693 ± 0.020 


0.095 ± 0.013 


1.355 ± 0.035 


0.160 ± 0.021 


0.49 


1.68 


0.001 


0.011 


0.001 


u 


- I 


1.624 ± 0.023 


0.114 ± 0.014 


2.505 ± 0.044 


0.196 ± 0.022 


0.49 


1.72 


0.001 


0.005 


0.001 


B 


- V 


0.707 ± 0.010 


0.046 ± 0.008 


0.904 ± 0.010 


0.045 ± 0.008 


0.51 


0.98 


0.001 


0.039 


0.001 


B 


- I 


1.642 ± 0.012 


0.058 ± 0.006 


2.056 ± 0.017 


0.079 ± 0.009 


0.49 


1.36 


0.001 


0.002 


0.001 


V 


- I 


0.931 ± 0.007 


0.038 ± 0.009 


1.148 ± 0.013 


0.052 ± 0.009 


0.49 


1.37 


0.001 


0.017 


0.001 
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Fig. 12. — Color distributions of GCs in NGC 1399, NGC 1404, and NGC 1387. Thick solid lines are smoothed histograms with Gaussian 
kernels of a{U — B) = 0.08, cr{B — V) = 0.04, and a-{V — I) = 0.04, respectively. The cr values are set to be comparable to the mean errors 
in each color. Filled gray histograms in the top panels show the color distributions of the brightest GCs (Vq < 21.5). Dotted histograms 
in the middle and bottom panels represent the color histograms of entire set of GCs within the extent of each galaxy, and solid histograms 
show the "corrected" distributions for the contamination by NGC 1399 GCs (see the text). 
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Fig. 13. — Variation in the red GC fractions (left) and tlie ratios of standard deviations between blue and red GC colors (right) for various 
colors. Black solid lines, red dashed lines, and blue dotted lines are for NGC 1399, NGC 1404, and NGC 1387, respectively. For NGC 1404 
and NGC 1387, the values were derived from the corrected color distributions for the contribution of NGC 1399 GCs in Figure [121 



are bluer at (U ~ B)o ~ —0.8 and foreground stars are 
at (f/ — B)o = —0.3 and 1.1 close to both ends of the 
GC range. The Vo-{B - V)o CMD (right panel) shows 
that many of the background galaxies have similar colors 
to the GCs, making it difficult to discriminate GCs from 
galaxies in this color. The foreground star sequences, 
particularly on the red side, are more clearly separated 
from GCs than in {U — B)q. This comparison corrobo- 
rates that the use of various color combinations including 
U -band colors provides a powerful tool for discriminating 
GCs from other objects. 

4.2. Color-Magnitude Diagrams 

It is interesting to note that the CMDs do not show 
bimodality for GCs brighter than Vq ~ 21.5. The 
disappearance was previously reported for NGC 1399 
(lOstrov et al.l Il998t iDirsch et alTl2003t iBlakeslee et"ar 



2012 ) , a few galaxies in th e Virgo Cluster (jMieske et al 
200a i Strader et al' '2006'), and o ther giant galaxies 



(jHarris et al. 2006; Faifer et al.l2011| ). The unimodal dis- 
tribution is often explained by merging of blue and red se- 
quences at the hig hest luminosities as a result of the blue 
tilt phenomenon ([Harris "eFaI|[2006[) . which is a trend 
for the blue GCs t o have redder color s at higher lumi- 
nositie s. However, lOstrov et al.l ()1998() and iForte et all 
(|2007| ) have reported that there is no blue tilt in the GC 
system of NGC 1399. Interestingly, our data reveal that 
the number of blue GCs rapidly drops at the brightest 
magnitudes (Vq < 21.5). This suggests that the unimodal 
distribution of brightest GCs is not the consequence of 
the blue tilt, but instead owing to the relative scarcity of 
blue, brightest GCs. 

4.3. Color Distributions 

Figure [T^ shows three representative color histograms 
of the entire GC system within the limiting radius in 
NGC 1399 (top), NGC 1404 (middle), and NGC 1387 



(bottom). The GCs fainter than the limiting magni- 
tudes and/or the ones located in the innermost region 
of each galaxy, where the completeness test is unreliable 
(see Section 2.3), are taken out from the subsequent anal- 
ysis. The thick solid lines are smoothed histograms with 
Gaussian kernels. The filled histograms in the top row 
are for the brightest GCs, as mentioned in the previ- 
ous section. For NGC 1399 (top row), GCs show bi- 
modal distributions for all colors considered in agree- 
ment with previous findings in C—Tl (J 3irsch_et al. 200^ 
IBassino et an[2f)06at Iftote et al.ll2007n ^r NGC 1404 
(middle) and NGC 1387 (bottom), the dotted histograms 
are for all GCs within the limiting radii of the galaxies. 
The solid histograms represent "corrected" distributions, 
for which the contribution of NGC 1399 GCs is statisti- 
cally subtracted from the original distributions and will 
be used in the further analysis. NGC 1404 GCs exhibit 
pronounced bimodality with an evident dip between two 
groups in all colors, consistent with the previou s stud - 
ies f or NGC 1404 ' s inne r GCs bv lGrillmair"eFan (|l999f ) 
and iLarsen et al.l (|2001| ). The clear separation between 
blue and red GCs for NGC 1387 GCs is also consis tent 
with what was found inC-Tlbyll assmo "ei~all (|2006bl. 
Blue GCs are much less abundant in NGC 1404 and NGC 
1387 than in NGC 1399. 

A closer scrutiny of Figure [12] reveals that the exact 
morphologies of the GC color distributions change de- 
pending on colors considered. The most distinctive case 
is given by the (t/-B)o distribution of NGC 1399. Com- 
pared to the {B — V)o and {V — /)o histograms, the 
(U — B)o distribution exhibits more prominent blue GC 
peak and have red GCs with weaker peak and larger dis- 
persion in color. The {U — B)q distribution shows spread 
in colors of red GCs that is nearly a factor of two larger 
than that of blue GCs, while the {B - V)o and {V - I)a 
histograms have similar variances for the red and blue 
GCs. The histograms for NGC 1404 and NGC 1387 also 
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show changes in their shape depending on colors. 

To be more quantitative in the bimodahty analysis, we 
used the Gaussian Mixture Modeling (GMM) code by 
iMuratov fc GnedinI (|2010( ). Table |6] presents the resuUs 
of the GMM analysis for six color combinations (([/ — 
B)o, (U - V)o, {U - I)o, {B - V)o, {B - /)o, and {V - 
/)o) for the three galaxies. It gives the color index, the 
mean (/j,) , and the standard deviation (ct) of blue and red 
GC colors, the red GC fraction (/,. = iV^d/A'^totai), and 
the ratio of the standard deviations between blue and 
red GCs [orloh). The last three columns summarize 
the probabilities of preferring a unimodal distribution 
over a bimodal distribution (p-values) derived based on 
the likelihood ratio test (LRT; x^): on the separation of 
the means relative to their variances (DD) and on the 
kurtosis of a distribution 

Figure [13] presents the results from the GMM analy- 
sis for all six color combinations. To quantify the shape 
of histograms, we use the red GC fractions (left) and 
the ratios of standard deviations between blue and red 
GC colors (right). NGC 1399 GC sample (black solid 
hues) docs not include GCs around NGC 1404 and NGC 
1387, and the GC sample of NGC 1404 (red dashed lines) 
and NGC 1387 (blue dotted) is corrected for NGC 1399 
GC contribution. For NGC 1399, we find that both the 
numbers and dispersions of blue and red GCs change 
significantly depending on colors that are used. The 
trend is more evident when comparing the ?7-band col- 
ors with the rest. For NGC 1404 and NGC 1387, the 
interesting features emerge after the correction for the 
contribution of NGC 1399 GCs. First, the red GC frac- 
tions in all colors remain constant regardless of the col- 
ors, as expected in Figure [12] showing distinct separa- 
tion between blue and red GCs. It is plausible that the 
early interactions of the galaxies with NGC 1399 have 
preferentially left the central, red GCs and later blue 
GCs have been accreted from outer region of NGC 1399. 
Second, the color dispersions vary significantly depend- 
ing on colors, following a similar pattern to the case of 
NGC 1399 for V - I, B - V, and B - I. For ;7-band 
colors, however, the patterns are different in the sense 
that they tend to be constant within the errors. The 
variation of dispersion as a function of colors is not ex- 
pected if the metallicity-to-color converting relations are 
straight. Instead, the slopes seem different between the 
blue and red parts of the metallicity-color relations and 
such an effe ct var i es syst ematically from color to color 
jYoon et aP [2001 1 20 11 alibi. ' 2012: C antiello fc Blakesled 
[2007t iChies-Santos et al.Tf20i2c iBlakeslee et al.ll20'T2[r 



4.4. Radial Variations of Color Bimodality 

Figure [13] displays the color distributions as a function 
of galactocentric radius for GCs in NGC 1399. The sub- 
sequent analysis is only valid for NGC 1399 because small 
number statistics prohibit us making any firm statement 
on the radial variations of GC colors for the two fainter 

^ Note that the kurtosis gives a necessary (but not suiHcicnt) con- 
dition of bimodality, and it can be used only as an additional check 
of other statistics ijMuratov fc GnedinI l2010l ) . For all galaxies, 
p(x^) from the LRT is < 0.01 in all colors, indicating the bimodal 
distributions are preferred over the unimodal ones. The p{DD) 
values lead to the same conclusion in all cases. The pikurt) values 
support the bimodal distributions in most cases, but favor a uni- 
modal distribution for {U-B)o of NGC 1399 with p{kurt) = 0.860. 



galaxies after the correction for the NGC 1399 GC con- 
tamination. In order to investigate the radial variation 
of color bimodality, we set up a series of radial bins con- 
taining approximately equal number 300) of GCs and 
perform the GMM analysis for each bin. We focus par- 
ticularly on three main parameters of the GMM analysis 
that characterize the bimodal distributions; the number 
ratio between blue and red GCs, the mean colors of the 
groups, and their color dispersions. 

Figure [15] shows the variations of color bimodality 
properties as a function of galactocentric radius. There 
is an obvious radial trend for the mean colors of both 
blue and red GCs getting bluer with increasing radius 
in all colors. The radial color gradient of red GCs is 
steeper than that of blue GCs. The red GC fraction in 
each color rapidly decreases with increasing galactocen- 
tric radius out to r 2± 4', and it remains fairly constant 
beyond the point. There is a weak radial trend in GC 
color dispersions, in that the ratio between dispersions 
of blue and red GCs increases with radius. 

5. SUMMARY AND CONCLUSIONS 

We have performed wide-field UBVI photometry of 
GCs in the central region of the Fornax cluster of 
galaxies with the Mosaic II CCD imager on the 4 m 
Blanco telescope at CTIO. This is one of the widest and 
deepest [/-band studies on extragalactic GC systems. 
The reduction was carried out with the DAOPHOT 
II/ALLFRAME package. Using two-color diagrams and 
the magnitude cut of My < —11, a total of 2037 GC 
candidates were selected among 12,134 objects detected 
in all four bands. Our estimate of the contamination by 
background galaxies and foreground stars in our sam- 
ple is in the range of 9% -16%. The sample is t/-band 
limited, and the completeness of our photometry is 50 % 
at Uq — 24.4 according to our artificial star tests. We 
provide the UBVI photometric catalog of the GC can- 
didates online. For the GC systems of NGC 1399, NGC 
1404, and NGC 1387, we have investigated the spatial 
distributions, CMDs, and color distributions. For NGC 
1399 GCs, the radial variations in the properties of bi- 
modality were examined as well. The main results are 
summarized as follows. 

1. We find an asymmetric distribution of GC candi- 
dates around NGC 1399. The overabundance of 
GCs at the southeast side of the inner region and 
toward the directions of NGC 1404 and NGC 1387 
in the outer region suggests that there were recent 
interactions of NGC 1399 with the galaxies and 
that NGC 1399 is not yet dynamically relaxed. 

2. We specify the radial extent of each GC system at 
the radius where the surface density profile begins 
to depart from the de Vaucouleurs fit. The limiting 
radii are ~3f4 for NGC 1404 and -3f7 for NGC 
1387. The GC system of NGC 1399 is spatially ex- 
tended beyond the field of view of our observation. 

3. We show that the GC systems in the three galaxies 
exhibit bimodal color distributions in all the colors 
considered in this study. For the brightest GCs in 
NGC 1399, the bimodality becomes weak or disap- 
pear depending on colors. 
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Fig. 14. — Histograms of {U — B)o, (B — V)o, and {V — I)o colors (left, middle, and right columns) for NGC 1399 GCs in various radial 
bins indicated in the upper-right corner within each panel. The radial bins are defined to contain approximately equal numbers (~ 300) of 
GCs. The histograms are expressed by two (i.e., blue and red) Gaussian normal distributions based on GMM analysis. The blue, rod, and 
black lines represent blue GCs, red GCs, and their sum, respectively. 
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Fig. 15. — Parameters of bimodality in {U — B)q, (B — V)q, and {V ~ I)o colors (left, middle, and right columns) as a function of 
galactocentric distance for NGC 1399 GCs. The top, middle, and bottom panels show the mean colors (fi), and /ir for blue and red GCs), 
the red GC fraction (/r), and the ratios of color dispersions between blue and red GCs (ar/tJb), respectively. 
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4. We find with GMM that the morphology of the 
color histograms changes significantly, depending 
on the colors used. The color dispersion of red GCs 
is larger than that of blue GCs in t/-band colors 
by more than a factor of two, while the red and 
blue GCs show similar dispersions in {B — V)o and 
(y — I)o- The number ratios between blue and red 
GCs vary depending on colors used as well. When 
corrected for the contribution of NGC 1399 GCs, 
the red GC fractions of NGC 1404 and NGC 1387 
remain constant regardless of the colors used. 

5. We suggest a possible scenario that could explain 
the presence of two discrete GC populations in 
NGC 1404 and NGC 1387: (1) the central, red 
GCs have been preferentially left through the early 
interactions of the galaxies with NGC 1399 and (2) 
later blue GCs have been accreted from outer re- 
gion of NGC 1399. 

6. We confirmed that the mean colors of blue and red 
GCs in NGC 1399 show a radial trend of becoming 



bluer with increasing galactocentric radius. The 
gradient of red GCs is steeper than that of blue 
GCs. 
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